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The foreseen capability to cover the far backward region at A Fixed-Target Experi¬ 
ment using the LHC beams allows to explore the dynamics of target fragmentation 
in hadronic collisions. In this report we briefly outline the required theoretical frame¬ 
work and discuss a number of studies of forward and backward particle production. 

By comparing this knowledge with the one accumulated in Deep Inelastic Scattering 
on target fragmentation, the basic concept of QCD factorisation could be investi¬ 
gated in detail. 
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I. INTRODUCTION 

In hadronic collisions a portion of the produced particle spectrum is characterised by 
hadrons carrying a sizeable fraction of the available centre-of-mass energy, the so-called 
leading particle effect. It is phenomenologically observed that for such hadrons their valence- 
parton composition is almost or totally conserved with respect to the one of inital-state 
hadrons [1]. In pp collisions, for example, protons, neutrons and lambdas show a significant 
leading particle effect. Moreover, for such semi-inclusive processes, the production cross 
section peaks at very small transverse momenta with respect to the collision axis, a regime 
where perturbative techniques can not be applied, giving insight on non-perturbative aspects 
of QCD dynamics in high energy collisions. 

Quite interestingly, the leading particle effect has been observed in Semi-Inclusive Deep 
Inelastic Scattering (SIDIS). At variance with the hadronic processes mentioned above, such 
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a process naturally involves a large momentum transfer. The presence of a hard scale enables 
the derivation of a dedicated factorisation theorem 2, 3] which ensures that QCD factori¬ 
sation holds for backward particle production in DIS. The relevant cross sections can then 
be factorised into perturbatively calculable short-distance cross sections and new distribu¬ 
tions, fracture functions, which simultaneously encode information both on the interacting 
parton and on the spectator fragmentation into the observed hadron. Despite of being 
non-perturbative in nature, their scale dependence can be calculated within perturbative 
QCD [ 4 ]. The factorisation theorem 2, 3] guarantees that fracture functions are universal 
distributions, at least in the context of SIDIS. 

Detailed experimental studies of hard diffraction at HERA have shown to support the 
hypothesis of QCD factorisation and evolution inherent the fracture functions formalism. 
Furthermore they led to a quite accurate knowledge of diffractive parton distributions Q- 
8], a special case of fracture functions in the very backward kinematic limit. For particles 
other than protons, proton-to-neutron fracture functions have been extracted from a pQCD 
analysis of forward neutron production in DIS in Ref. { 9 ]. A set of proton-to-lambda fracture 
functions has been obtained by performing a combined pQCD fit to a variety of semi-inclusive 
DIS lambda production data in Ref. [10]. 

As theoretically anticipated in Ref. 2,[ll|, 12] and experimentally observed in hard diffrac¬ 
tion in pp collisions at Tevatron jl3, 14], QCD factorisation is violated for fracture functions 
in hadronic collisions. On general grounds, it might be expected, in fact, that the dynamics 
of target-remnants hadronisation is affected by the coloured environment resulting from the 
scattering in a rather different way with respect to the Deep Inelastic Scattering case. 

Nonetheless, the tools mentioned above allow to investigate quantitatively particle pro¬ 
duction mechanisms in the very backward and forward regions, to test the concept of fac¬ 
torisation at the heart of QCD and to study the dependencies of factorisation breaking upon 
the species and the kinematics of the selected final state particle. 

This physics program could be successfully carried on at A Fixed-Target Experiment 
using the LHC beams jl5]. Novel experimental techniques are, in fact, available to extract 
beam-halo protons or heavy-ions from LHC beams without affecting LHC performances. 
Such a resulting beam would be then impinged on a high-density and/or long-lenght fixed 
target, guaranteeing high luminosities. Furthermore and most importantly for the physics 
program to be discussed here, the entire backward hemisphere (in the centre-of-mass system 
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of the collision) would be accessible with standard experimental techniques allowing high 
precision studies of target fragmentation. Althought measurements of particle production 
in the very forward region (close to the beam axis) might be challenging experimentally due 
to the high particle densities and large energy flow, the installation of dedicated detectors, 
like forward neutron calorimeters and/or proton taggers, could further broaden the physics 
program oulined above giving access to the beam fragmentation region. 

The paper is organised as follows. In Sec. Owe first give a brief theoretical introduction 
on the fracture functions formalism and to higher order corrections to the semi-inclusive 
Drell-Yan process. In Sec. II I II we outline different analyses which could performed at AF- 
TER@LHC with special focus on single hard diffraction. In Sec. IIVI we summarise our 
results. 


II. COLLINEAR FACTORISATION FORMULA 

Fracture functions, originally introduced in DIS, do depend on a large momentum trans¬ 
fer. Therefore, in order to use them in hadronic collisions, a hard process must be selected. 
We consider here the semi-inclusive version of the Drell-Yan process, 

tfr(Pi) + H 2 (P 2 ) H(h ) + 7 *(g) + X, (1) 

in which one hadron H is measured in the final state together with a Drell-Yan pair. In 
such a process the high invariant mass of the lepton pair, q 2 = Q 2 , allows the applicability 
of perturbative QCD, while the detected hadron H can be used, without any phase space 
restriction, as a local probe to investigate particle production mechanisms. 

The associated production of a particle and a Drell-Yan pair in term of partonic de¬ 
grees of freedom starts at 0(a s ). One of the contributing diagrams is depicted in Fig. (JT|) . 
Assuming that the hadronic cross-sections admit a factorisation in term of long distance 
non-perturbative distributions and short distance perturbative calculable matrix elements 
for the partonic process i(pi) + —> l(k ) + 7 *(<?), predictions based on perturbative 

QCD are obtained convoluting the relevant partonic sub-process cross-sections, dcW - ^ 7 *, 
with parton distribution functions, /)/ h x and fj/H 2 i and fragmentation functions, D H ^ 1 . The 
hadronic cross section, at centre of mass energy squared s = (P\ + P 2 ) 2 , can be symbolically 
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FIG. 1: Example of diagram contributing to hadron production in the central fragmentation region 
to order 0(a s ) in eq. (El). 


written as 



da H ’ C ’W 

dQ 2 dz 



1 ) f? ] (.x 2 ) D H/l (z/p ) 
P 


da ij ^* 
dQ 2 dp ’ 


( 2 ) 


where the convolutions are over the momentum fractions of the incoming and outgoing 
partons. The partonic indeces i , j and l in the sum run on the available partonic subpro¬ 
cesses. The superscripts label the incoming hadrons and the presence of crossed terms is 
understood. This type of factorised hadronic cross section is expected to hold for hadrons 
produced at sufficiently high transverse momentum and it is widely and successfully used to 
compute cross sections for large momentum transfer processes in hadronic collisions. The 
Lorentz-invariant variable z in eq. ([2]) is defined by 


z = 


2 h ■ (Pi + P 2 ) _ 2 


(3) 


s \/s 

In the hadronic centre-of-mass frame, where the second identity holds, z is just the observed 
hadron energy, E* H) scaled down by the beam energy y^/2. The variable p, appearing in 
eq. ([2]), is its partonic equivalent. Within this production mechanism, the observed hadron 
H is generated by the fragmentation of the final state parton /, and for this reason we ad¬ 


dress it as central. The amplitudes squared 


18], however, are singular when the transverse 


momentum of the final state parton vanishes. In such configurations, the parent parton 
l of the observed hadron H is collinear either to the incoming parton i or j. As these 
phase space region are approached, perturbation theory looses its predictivity. This class of 















5 



FIG. 2: Pictorial representation of the parton model formula, eq. m, for the associated production 
of a Drell-Yan pair and a particle in the target fragmentation regions. 

collinear singularities escape the usual renormalisation procedure which amounts to reab¬ 
sorb collincar divergences into a redefinition of bare parton distribution and fragmentation 
functions. Such singularities are likely to appear in every fixed order calculation in the same 


kinematical limits spoiling the convergence of the perturbative series. In Refs 


16l |17j a gen¬ 


eralised procedure for the factorisation of such additional collincar singularities is proposed. 


The latter is the same as the one proposed in Deep Inelastic Scattering 


19] where the same 


singularities pattern is also found, confirming the universality of collinear radiation between 
different hard processes. Such a generalised collinear factorisation makes use of fracture 
functions. These distributions obey DGLAP-type evolution equations which contain an 
additional inhomogeneous term resulting from the subtraction of collinear singularities in 


the target-fragmentation region 


19]. Such equations allow to resum the corresponding 


large logarithmic corrections to all orders in perturbation theory. Bare fracture functions, 
M i 7 (x,z), describe the hadronization of the spectators system in hadron-induced reac¬ 

tions. They express the conditional probability to find a parton i entering the hard scattering 
while an hadron H is produced with fractional momentum z in the target fragmentation 
region of the incoming hadron H\. 

The use of fracture functions allows particle production already to O(a 0 s ), since the 
hadron H can be non-pertubatively produced by a fracture function itself. Therefore the 
lowest order parton model formula can be symbolically written as 


dcr H ’ T,( '°> f dx i f dx2 r . r fii , , P yi\, \ *Ai\, N nil, , n dd 1 ^ 1 * . 

/ ~ / — [ M i (si, z ) f\ (afr) + M \ (g 2 , z) f) (si)] (4) 

i,j 


dQ 2 dz 
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FIG. 3: Example of diagram contributing to 0(a s ) corrections in the target fragmentation region, 
eq.$5j). 


and it is sketched in Fig. ((2]). The superscripts in eq. flj]) indicate from which incoming 
hadron, Hi or H 2l the outgoing hadron H is produced through a fracture functions. In order 
to complete the calculation to 0(a s ) accuracy we should consider higher order corrections 
to eq. »• Since in this case the hadron H is already produced by fracture functions, final 
state parton radiation should be integrated over and the resulting contribution added to 
virtual corrections. One of the contributing diagrams is depicted in Fig. (J3]). The general 
structure of these terms is 


da H ' T ' ID 
dQ 2 dz 


oc 


E 




— I — [ M l 1] ( x h z ) ff\ x 2 ) + M f\ x 2 , z) ff\x 1 )] —— . (5) 

X\ J X2 CiKcJ 


We refer to them as to the target fragmentation contributions. Their calculation is, a part 
from minor differences in kinematics, completely analogous to the one of the inclusive Drcll- 


Yan case. The factorisation procedure, first elaborated in Ref. Hi in the context of SIDIS, 
amounts to substitute in eq. (HI) the bare fracture and parton distributions functions with 

12]. Renormalised parton distributions and fracture functions 


their renormalised version 


16 


homogeneous terms do cancel, as in the inclusive Drcll-Yan case, all singularities present in 
eq. (|5|) . The additional singularities in eq. ([2]) are cancelled by the combination of parton 
distributions and fracture functions inhomogeneous renormalisation terms. Adding all the 
various contributions, the resulting hadron-/y integrated cross section, up to order 0(a s ), is 
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then infrared finite 



and can be simbolically written as 


da 11 do 
dQ 2 dz N c s 


EN. 


[i] 


if + Mf 1 


/j 11 ]( 1 + s^)+ 




2tt 

(Tq Ol s 


N r s 27r 


- V/ 


[ 1 ] 


V 

J 3 


[ 2 ] 


D H / l ®K\ j , (6) 


mV 


where cr 0 = 4:7ral m /3Q 2 and N c is the number of colors. We refer to the previous equation 
as to the collincar factorisation formula for the process under study. The next-to-leading 
order coefficients C y ' and have been calculated in Ref. jl7|, making the whole calculation 
ready for numerical implementation. 

We stress, however, that our ability to consistently subtract collincar singularities in such 
a semi-inclusive process is a necessary but not sufficient condition for factorisation to hold 
in hadronic collisions. The one-loop calculation outlined above in fact does involve only 
the so-called active partons. It completely ignores multiple soft parton exchanges between 
active and spectators partons, whose effects should be accounted for in any proof of QCD 
factorisation. Therefore there is no guarantee that fracture functions extracted from SIDIS 
can be successfully used to describe forward or backward particle production in hadronic 
collisions. Reversing the argument, such a comparison may instead offer new insights on 
non-perturbative aspects of QCD and to the breaking of factorisation. 


III. SINGLE HARD DIFFRACTION AT AFTER@LHC 

As an application of the formalism presented in the previous sections we will consider 
single hard diffractive production of a Drell-Yan pair 

Pi(Pi) + P 2 (A) p(P) + 7*(<z) + *, (7) 

where we have indicated in parathesis the four momenta of the relevant particles. We 
present in the following cross sections differential in the virtual photon variables. The 
subsequent decay of the virtual photon into a lepton pair can be easily included so that 
realistic cuts on leptons rapidity and transverse momentum can then be applied. We consider 
the AFTER@LHC kinematic setting in which a 7 TeV proton beam collides on a fixed target 
proton leading to a centre-of-mass energy of yfs = 115 GeV. We consider the projectile 
proton pi moving in the positive z direction and p 2 at rest in the laboratory. The diffractively 








produced proton p has in general almost the incoming projectile proton p\ energy and very 
small transverse momentum as measured with respect to the collision axis. The detection 
of such fast protons will in general require the installation of forward proton taggers. The 
lepton pair instead will be measured by the main AFTER@LHC detector. This kinematical 
configurations is pictorially represented in right plot of Fig. (J2j). 

Diffractive processes has been intensively analysed in DIS at HERA ep collider, revealing 
their leading twist nature. From scaling violations of the diffractive structure functions 


and dijet production in the final state [6, 7] quite precise diffractive parton distributions func¬ 
tions (dPDFs) have been extracted from HERA data, which parametrise the parton content 
of the color singlet exchanged in the t-channel. The comparison of QCD predictions for sin¬ 
gle diffractive hard processes based on diffractive parton distributions measured at HERA 
against data measured at Tevatron h], 14] (i/s = 1.96 TeV), adopting a factorised ansatz 
as in eq. ([3]), has indeed revealed that these processes are, not unexpectedly ll|, |l2], signif¬ 
icantly suppressed in hadronic collisions. This conclusion persists even after the inclusion 


of higher order QCD corrections 


20 ]. Complementing these results with the forthcoming 


ones at LHC at higher at yfs = 13 TeV and the ones from AFTER@LHC at yfs = 115 
GeV will give insight on the energy dependence of the socalled rapidity gap survival (RGS) 
probability in a wide range in yfs. Since the theoretical computation of the RGS factor is 
highly model dependent, we decided not to include it in our predicted cross sections which 
must be considered therefore as upper bounds. 

Diffractive parton distributions ff are in general proton-to-proton fracture functions Mi. 
They depend upon the final state proton fractional energy loss, Xip = 1 — z with z given 
in eq. ([3]), the fractional momentum of the interacting parton with respect to the pomeron 
momentum, (3 = x/xjp and the virtuality Q 2 . In general fracture functions may depend 
also upon the invariant momentum transfer t — (P — Pi) 2 at the proton vertex [21]. In all 
diffractive structure functions measuraments at HERA, out of which dPDFs are determined, 
t is integrated over up to some t max <C Q 2 . In this case dPDFs obey ordinary DGLAP 
evolution equations |22| as their extended, f-dependent, version [21]. In the present paper 


we use dPDFs form Ref. [5[ which are defined by |t| < 1 GeV 2 . Since they are extracted 
from large rapidity gap data where the proton is not directly measured, they contain a 
contribution (23%) from the socalled proton dissociation contribution. In order to use 
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FIG. 4: Left: Double differential cross sections for the production of a Drell-Yan pair of mass 
Q 2 = 100 GeV 2 . Blue error bands represent theoretical errors estimation, as described in the text. 
Right: Double differential cross sections times Q 4 for three different invariant masses. 


dPDFs in the present context we first note that 


Mi(x i, z,Q )= x IP /,; (0, x IPl Q 


( 8 ) 


The extra factor Xjp comes from the jacobian of the change X\ = fdxip. We then rearrange 
eq. (HJ) in terms of new variables obtaining 

dcr D do rl 

dQ 2 dxip N c s j T j q 


h, J Y x, P , A) f^, A) . 


( 9 ) 


with r = Q 2 /s. For simplicity we consider here leading order formulas but the extension to 
higher order is straightforward. In eq. (j^J) we use parton distribution functions from Ref. 23]. 
We show explicitely the dependence of fracture and parton distributions functions upon the 
factorisation scale, p 2 F . Predictions are obtained with this scale set to p 2 F = Q 2 . Theoretical 
errors associated with higher order corrections are instead estimated varying such scale in 
the range p 2 F = 1/4 Q 2 and p 2 F = 4 Q 2 . 

In Fig. (j3J) we present predictions for the xjp distribution. In left plot we consider a 
Drell-Yan pair of mass Q 2 = 100 GeV 2 . The distribution shrinks as lower x/p-values are 
approached whereas, from hard diffraction at HERA, it is well known that diffractive cross 
sections rise as an inverse power of xj P . Such effect therefore is then attributed to phase 
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FIG. 5: Left: Double differential cross sections for the production of a Drell-Yan pair at xjp = 0.05. 
Blue error bands represent theoretical errors estimation, as described in the text. Right: Double 
differential cross sections times Q 4 for three different xip-values. 

space threshold effects. The Drell-Yan invariant mass constraint can be rewritten in the 
diffractive case as Q 2 = /3xjpx 2 s, which can be cast (for /3 —* 1 and x 2 —> 1) in upper bound 
on the invariant mass Q 2 < Xjps at fixed Xip and s. This hypothesis is further supported in 
the right plot of Fig. ()4]), where differential distributions are presented for three values of Q 2 . 
The former is multiplied by Q 4 to compensate the fast fall off of the electromagnetic cross 
section. The lowest values of xip are then accessed only by lowering the invariant mass of 
the pair. We note that, even considering the maximum value of Xjp = 0.1, single diffractive 
production of and Z is beyond the kinematic reach at AFTER@LHC. In the left panel 
of Fig. (J5]) we present the prediction for the Q 2 distribution at a fixed value of xip = 0.05. 
The cross section, as expected, is fast falling as an inverse power of Q 2 . The Q 2 distribution 
is particularly instructive since it allows to study the possible dependence of the RGS factor 
on Q 2 and therefore to determine the underlying dynamics. In the right panel of Fig. ([5]) we 
present the (^-differential cross section again multiplied by the factor Q 4 . In this way all the 
Q 2 dependence is accounted for by that of fracture and parton distributions. These curves 
and the corresponding slopes, however, can not be readily interpreted as genuine results of 
QCD evolution of fracture and parton distributions functions because of the threshold effect 
mentioned above appearing at such moderate values of i/s. 
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FIG. 6: Left: Triple differential cross sections for the production of a Drell-Yan pair at xjp = 0.05 
and of mass Q 2 = 100 GeV 2 . Blue error bands represent theoretical errors estimation, as described 
in the text. Right: Triple differential cross sections for three different xjp-values. 


By changing variable from (3 to the virtual photon centre-of-mass rapidity, y 

d T n ,CM 


CM 


P = 


Xip 


/— 

, X 2 = V r < 


-y 


( 10 ) 


eq. 0 can be further manipulated to give the three-differential cross section 

dcr D 


<?o 


dQ 2 dx IP dy CM N c s 


^2 e l X Ipf°(Pi XlP ’ Vf) fq( X 2, Tf) 


( 11 ) 


q,q 


The rapidity range for diffractive Drell-Yan production reads 


In y/r < y CM < In (12) 

xip 

which, as expected, turns out to be asymmetric given the kinematic constraint X\ < Xip. The 
rapidity range for the inclusive Drell-Yan case is recovered simply setting Xip = 1. The ra¬ 
pidity distributions is particularly sensitive to the shape the diffractive parton distributions. 
This distribution will be useful to investigate any possible kinematic dependence of the RGS 
factor. In the left panel of Fig. ([6]) we present the centre-of-mass rapidity distribution at 
fixed Q 2 = 100 GeV 2 and xip = 0.05. In this frame the distribution is shifted at negative 
values of y CM . Therefore, on average, the parton originating from the target proton carries 
more momentum than the one originating from the pomeron. Since the rapidity is additive 
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FIG. 7: Left: Rapidity distributions for inclusive and diffractive Drell-Yan of mass Q 2 = 100 GeV 2 . 
Blue error bands represent theoretical errors estimation, as described in the text. Right: diffractive 
to inclusive Drell-Yan rapidity distributions ratio. 

under boost along the collision axis we may easily boost the y CM to the laboratory frame 
by using 

/■“ = + y CM with 0=^1-!^ (13 ) 

with m p the proton mass. In the AFTER@LHC kinematics this implies a rapidity shift 
Ay = y Lab — y CAI = 4.8. The rapidity distributions in the laboratory frame for a Drell- 
Yan pair of mass Q 2 = 100 GeV 2 and for three different xjp values are presented in the 
right panel of Fig. (J6]). One may notice from the plot that for increasing Xjp the Drell- 
Yan pair spans a wider rapidity range and the corresponding spectrum is increasingly more 
forward. It might be useful to discuss the single diffractive Drell-Yan pair production in 
conjuction with the analogous inclusive process. Such a comparison is presented for the 
centre-of-mass rapidity distributions in the left plot of Fig. (J7J) for a common Drell-Yan 
pair of mass squared Q 2 = 100 GeV 2 and, for the diffractive case, integrated in the range 
10~ 3 < xip < 10” 1 . The rapidity distributions in the single diffractive case is strongly 
asymmetric whereas in the inclusive case it is symmetric around y CAl = 0. This effect is 
primarily due to the different kinematics of the two processes and to the different fractional 
momentum distributions between parton and fracture distributions. In the right plot of 
Fig. (HD the ratio between the two previous distributions is presented. The ratio gives 
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direct information on the suppression factor between the single diffractive to the inclusive 
process, assuming a factorised expression for the former, eq. ( 1111 ) . Such ratio might be 
convenient from the experimental side since many lepton detection systematics will cancel. 
On the theoretical one it is expected to be more stable against the inclusion of higher order 
corrections. In the present case, for example, the factorisation scale is simultaneously varied 
both on the numerator and denominator resulting in a reduced theoretical error band with 
respect to the one obtained for absolute cross sections. 

We wish to end this section we a brief overview of other possible applications of the 
proposed formalism. A completely analogous program can be performed for the associated 
production of forward neutron and a Drcll-Yan pair, p + p —» n + 7 * + A". The production of 
forward neutron in DIS at HERA has shown a leading twist nature. From scaling violations 
of the semi-inclusive neutron structure functions a set of proton-to-neutron fracture functions 
set has been extracted from data in Ref. jsi] which can be used to predict forward neutron 
rate in hadronic collisions. As in the case of hard diffraction, both physics programs would 
highly benefit from the installation of a dedicated instrumentation for the measuraments 
of fast neutrons and protons quite close to the beam axis. Measuraments in the forward 
region, althought problematic experimentally, give in fact direct access to the study of the 
beam fragmentation region. 

As a third application we consider hyperon production associated with a Drcll-Yan pair, 
p + p —> V + 7 * + A", where V generically indicates either a A 0 or A 0 hyperon. At very 
low transverse momentum, A 0 longitudinal momentum spectrum should show a significant 
leading particle effect, which can be predicted, assuming factorisation as in da H,T , by the 
proton-to-lambda fracture functions set obtained from a fit to SIDIS data in Ref. 10]. On 
the other hand A 0 spectrum in the same kinematical conditions should instead show almost 
no leading particle effect, giving access to the proton-to-A 0 fracture functions. We note, 
in general, that the particle-to-antiparticle fracture function is indeed an interesting and 
almost unknown distribution. On the other hand, if one considers A 0 or A 0 at sufficiently 
large transverse momentum, their combined analysis, described by da H,c , should allow the 
investigation of parton hadronisation into hyperons in the QCD vacuum as parametrised by 
fragmentation functions. 

As a last application we consider the associated production of one particle and a Drell-Yan 
pair in the context of multi-parton interactions. The latter process has already been used to 
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investigate the contamination of the so-called underlying event 
been successfully used to study underlying event properties 


24] to jet observable and has 
25] . If the detected hadron is 


measured at sufficiently large transverse momentum, the latter constitutes a natural infrared 
regulator for the partonic matrix elements. In this kinematics conditions we also expect a 
rather small contributions from fracture functions. Therefore the central term, dcr 11,0 ’^, 
can be used to estimate the single parton scattering contribution to the process. The latter 
might be considered as the baseline to study the contributions of double (or multiple) parton 
scattering contributions to the same final state, where, for example, the primary scatter 
produces a Drcll-Yan pair while the secondary one produces the detected hadron H. 


IV. CONCLUSIONS 

We have briefly reviewed a perturbative approach to single particle production associated 
with a Drell-Yan pair in hadronic collisions. On the theoretical side we have shown that the 
introduction of fracture functions allows a consistent factorisation of new class of collinear 
singularities arising in this type of processes. The factorisation procedure coincides with the 
one used in DIS confirming, as expected, the universal structure of collinear singularities 
and supports the proposed collinear factorisation formula. On the phenomenological side 
we have outlined some areas in which the formalism can be fully tested. In particular, 
focusing on the AFTER@LHC kinemtic range, we have discussed in some detail the single 
diffractive production of virtual photons. The study of such a process might improve our 
understanding of non-perturbative aspects of QCD and it allows to explore in detail the 
nature of factorisation breaking at intermediate energies. 
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